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Spread Spectrum Signal Processing 

5 THE BACKGROUND OF THE INVENTION AND PRIOR ART 

The present invention relates generally to processing of spread 
spectrum signals by means of vector based algorithms. More 
particularly the invention relates to a method of receiving spread 
spectrum signals according to the preamble of claim 1 and a 
10 signal receiver according to the preamble of claim 28. The 
invention also relates to a computer program according to claim 
26 and a computer readable medium according to claim 27. 

Spread spectrum technology is becoming increasingly important 
in many areas, of which cellular communication systems and 

15 global navrgatlon satellite systems (QNSS) represent two 
important examples. Moreover, the large variety of transmission 
standards creates a demand for hybrid- or general-purpose 
receivers. For instance, both the cellular standards cdma2000 
and WCDMA involve transmission of spread spectrum signals. 

20 However, due to differences In the air interfaces one and the 
same terminal cannot be used in the two systems. Instead, a 
dedicated terminal must be employed in each system. Alter- 
natively, a dual mode terminal must be designed, which Includes 
two separate transceiver chains. Correspondingly, a navigation 

25 receiver adapted for one GNSS, say the Global Positioning 
System (GPS; U.S. Government), is not able to receive signals 
from a satellite that belongs to a different GNSS. such as the 
Galileo system (the European programme for global navigation 
services) or the Global Orbiting Navigation Satellite System 

30 (GLONASS; Russian Federation Ministry of Defense). 
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In order to make this possible, a multi-mode receiver must be 
designed. However, including multiple receiver chains in a 
single device is not only expensive, it also renders the unit bulky 
and heavy, particularly if more than two signal formats are to be 
.5 processed. Therefore, a programmable software receiver 
solution is desired whose signal processing principles may be 
altered according to which signals that shall be received and 
processed. 

Various software solutions are already known for processing 
10 GNSS signals. The patent document WO02/50561 is one 
example in which a GPS receiver tracking system for guiding 
missiles is described. The receiving station here includes one or 
more processing stages for receiving and processing GPS data. 
A control signal controls whether a quantity of processing stages 
IS should be increased or decreased in response to the result of a 
correlation operation. The tracking and recovery of timing in- 
formation is entirely performed in software. However, the pro- 
cessing involves reaMime calculation of relatively complex 
Fourier transforms, and therefore requires considerable pro- 
20 cessing resources. Thus, a receiver of this type is not parti- 
cularly suitable for small sized units, such cellular phones or 
portable GNSS-receivers. 

The article "TUPM 12.4: Software Solution of GPS Baseband 
Processing". ULSI Laboratory, Mitsubishi Electric Corporation, 

25 IEEE, 1998 by Asai T.. et al. outlines the principles for a soft- 
ware GPS receiver using an embedded microprocessor, which 
enables signals from up to eight satellites to be demodulated at 
around 40 MIPS (millions of / or Mega instructions Per Second). 
A 2 MHz GPS baseband signal is here fed to the software radio 

30 receiver, which performs spectrum despreading and mixing; 
synchronous demodulation, satellite acquisition and tracking; as 
well as phase compensation. An algorithm is proposed through 
which the intermediate frequency (IF) signal Is 1-bit sampled 
and where downconversion is performed before despreading. 

35 Resulting data in the form of 32 samples are continuously 
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processed in parallel. A DRAM (Dynamic Random Access 
Memory) stores a particular C/A (Coarse Acquisition) code table 
for each satellite from which signals are received. All other data 
and Instructions are stored in a cache memory in order to 
5 accomplish a high-speed processing. Although this strategy is 
declared to result in a reasonably efficient impiementatton, the 
document lacks a specific description as to how the dispreading 
and the correlation should be effected in order to, in fact, obtain 
this effect. 

10 Akos. D. et al.. "Tuning in to GPS - Real-Time Software Radio 
Architectures for GPS Receivers", GPS World, July 2001 
describes a receiver architecture througli which IF signal 
samples are fed directly from a radio front-end to a program- 
mable processor for continued processing. The article mentions 

15 the possibility of using single instruction multiple data (SIMD) 
instructions to process multiple data samples in parallel. How- 
ever, portability- and flexibility problems are recognized, and 
again, there is no teaching as to how an efficient parallel pro- 
cessing could actually be accomplished. 

20 Dovis, F. et ai., "Design and Test-Bed Implementation of a 
Reconfigurabie Receiver for Navigation Applications", Electro- 
nics Department, Ppiitecnico di Torino. Navigation Signal Ana- 
lysis and Simulation Group, Spring of 2002 relates to the design 
of a reconfigurabie GNSS receiver which is capable of fusing 

25 data from two or more different GNS8:s. The document sketches 
an architecture which, in addition to a radio front-end, includes a 
Field Programmable Gate Array (FPGA) and a Digital Signal 
Processor (DSP). The authors address various computational 
load issues. However, they do not present an algorithm which 

30 fulfills the identified requirements. 

Hence, the prior art Includes a number of examples of software- 
based GNSS-receivers. Nevertheless, there is yet no distinct 
teaching of a highly efficient solution which is suitable for imple- 
mentation in software and has capabiiities that are at least in par 



with those of today's AStC-based solutions for receiving and pro- 
cessing spread spectrum signals in real-time (ASIC = Application 
Specific integrated Circuit). 



SUMMARY OF THE INVENTION 

5 The object of the present invention is therefore to provide a 
solution for receiving and processing spread spectrum signals, 
v^ich solves the problems above and thus to present a truly 
software adapted strategy that may be implemented in a 
physically small and power efficient device. 

10 According to one aspect of the invention the object is achieved 
by a method of receiving spread spectrum signals as initially 
described, which is characterized by comprising a preparation 
for the correlation step, wherein the preparation takes place 
before the receiving the continuous signal. This preparation 

15 involves pre-generating a multitude of code vectors, which each 
represents a particular code sequence of the at least one signal 
source specific code sequence. Moreover, according to the 
invention, the correlation step involves multiplying at ieast each 
vector In a sub-group of the code vectors with at least one 

20 vector, which is derived from the data word. 

This strategy is advantageous because the proposed vector 
approach makes it possible for a digital processor (e.g. a micro- 
. processor) to process multiple signal samples in parallel during 
each clock cycle and thus make very efficient use of the 
25 processor. Furthermore, by pre-generating the code vectors 
valuable processing capacity is saved, which in turn renders It 
possible to deal with a sufficient amount of input data per time 
unit in order to. for example track navigation signals that are 
transmitted by the satellites in a navigation satellite system. 

30 According to a preferred embodiment of this aspisct of the 
invention, each code vector represents a particular signal 
source specific code sequence, which is sampled at the basic 
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sampling rate and quantised with the quantising process that is 
used to produce the level-discrete sample values. Preferably, 
the code vectors are also further adapted to the format of the 
incoming data signal in a way which renders the following 
5 processing efficient. 

According to another preferred embodiment of this aspect of the 
invention, the at least one signal source specific code sequence 
represents so-called pseudo random noise. This type of code 
sequences is beneficial in environments where two or more 

10 different signals are modulated onto the same carrier frequency. 
Namely, the different pseudo random noise signals are ortho- 
gonal (or at least almost orthogonal) to each other. Conse- 
quently, a first information signal spread by means of a first 
. pseudo random noise signal causes a minimal deterioration of a 

IS second information signal which is spread by means of a second 
pseudo random noise signal, and vice versa. 

According to another preferred embodiment of this aspect of the 
invention, the receiving step involves down conversion of an 
incoming high-frequency signal to an intermediate frequency 

20 signal. The high-frequency signal is presumed to have a spec- 
trum which is symmetric around a first frequency and the inter- 
mediate frequency signal is presumed to have a spectrum which 
is symmetric around a second frequency, which is considerably 
. lower than the first frequency. Hence, the receiving step 

25 transforms the information carrying signal to a baseband which 
is technically more easy to handle than the band at which the 
original incoming signal is located. 

According to another preferred embodiment of this aspect of the 
invention, the method includes the following steps. First, a 
30 maximum frequency variation of the second frequency due to 
Ooppler-effects is determined. Then, a Doppler frequency interval 
around the second frequency is defined. The Doppler frequency 
interval has a lowest frequency limit equal to the difference 
between the second frequency and the maximum frequency 
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variation, and a highest frequency limit equal to the sum of the 
second frequency and the maxirnum frequency variation. After that, 
the Doppler frequency interval is divided into an integer number of 
equidistant frequency steps, and subsequently a frequency 
5 candidate vector is defined for each frequency step. This way of 
dividing the spectrum is advantageous because it provides a 
flexible modeling of the intermediate carrier frequency and its 
variations. 

According to another preferred embodiment of this aspect of the 
10 invention, an integer number of initial phase positions is 
determined for the intermediate frequency signal. The integer 
number thus represents a phase resolution, such that a 
relatively high number corresponds to a comparatively high 
phase resolution, whereas a relatively low number corresponds 
15 to a comparatively low phase resolution. Furthermore, a carrier 
frequency>phase candidate vector is defmed for each combi- 
nation of carrier frequency candidate vector and initial phase 
position. The different frequency-phase candidate vectors are 
beneficial, since they allow the determination of a very accurate 
20 estimate of frequency* and phase characteristics of a received 
signal, which In turn vouches for a demodulated signal having a 
high quality. 

According to yet another preferred embodiment of this aspect of 
the invention, tiie number of elements in each carrier frequency- 

25 phase candidate vector is determined. Subsequently, the carrier 
frequency-phase candidate vectors are stored according to a 
data format, which is adapted to performing correlation 
operations between the at least one vector derived from the data 
word and a segment of a carrier frequency-phase candidate 

30 vector. This adaptation Is advantageous because thereby 
several signal samples may be processed In parallel by means 
of comparatively simple operators. Preferably, the adapting of 
the data format involves adding at least one element to each 
segment of the carrier frequency- phase candidate vector, such 

35 that the segment attains a number of elements which Is equal to 
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the number of elements in each of the at least one vector that is 
derived from the data word. Hence, it Is namely possible to 
process a segment of the carrier frequency-phase candidate 
vector together with one of the at least one vector by either a 
5 so-called Single Instruction Multiple Data (SIMD)-operation or a 
logical XOR-operation. 

According to stilt another preferred embodiment of this aspect of 
the invention, the method also includes the following steps. 
First, a maximum variation of the code rate due to Doppler- 

10 effects Is determined. Then, a Doppler rate interval Is defined 
around a center code rate. The Doppler frequency interval has a 
lowest code rate limit equal to the difference between the center 
code rate and the maximum code rate variation, and a highest 
frequency limit equal to the sum of the center code rate and the 

15 maximum code rate variation. In analogy with the intermediate 
frequency spectrum, the Doppler rate interval is divided into an 
Integer number of equidistant code rate steps, and a code rate 
candidate is defined for each codel rate step, in further analogy 
with the intermediate frequency, an integer number of possible 

20 initial code , phase positions is determined for each code rate 
candidate. The integer number thus represents a code phase 
resolution, such that a relatively high number corresponds to a 
comparatively high code phase resolution, whereas a relatively 
low number corresponds to a comparatively low code phase 

25 resolution. Typically, the spread spectrum code modulation is 
less sensitive to distortion caused by Doppler effects than the 
. frequency modulation. Therefore, the code rate steps may be 
relatively large. However, a certain degree of Doppler shift esti- 
mation is desirable to attain a demodulated signal of high quality. 

30 According to yet another preferred embodiment of this aspect of 
the invention, a set of code vectors is generated for each signal 
source specific code sequence by sampling each code rate- 
phase candidate vector with the basic sampling rate, and thus 
produce a corresponding code vector. Again, this simplifies a 

35 parallel processing of several signal samples by means of 
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comparatively simple operators. 

According to still another preferred embodiment of this aspect of 
the invention, a modified code vector is generated on basis of 
each code vector by: copying a particular number of elements 

5 from the end of an original code vector to the beginning of the 
modified code vector, and copying the particular number of 
elements from the beginning of the original code vector to the 
end of the code vector. This extension of the code vector is 
highly advantageous because It allows extraction of tracking 

10 parameters through early-, prompt- and late techniques by 
performing a simple tranislation of local copies of the signal 
source specific code sequence. 

According to still another prefenred embodiment of this aspect of 
the invention, a set of modified code vectors is stored for each 

15 signal source specific code sequence. Each modified code 
vector here contains a number of elements, which represent a 
sampled version of at least one full code sequence. As 
mentioned above, the modified code vector also includes a 
repetition of the beginning and the end of the sequence. A 

20 particular modified code vector is defined for each combination 
of code rate candidate and code phase position. These pre- 
generated modified code vectors save considerable processing 
capacity in the processor, since they need not be calculated in 
real time, but may simply be acquired from a memory. 

25 According to yet another preferred embodiment of this aspect of 
the invention, the data format of the modified code vectors Is 
adapted with respect to the data format of the at least one 
vector that has been derived from the data word, such that a 
modified code vector and one of the at least one vector may be 

30 processed jointiy by a SlMD-operation or an XOR-operation. 
Naturally, this Is advantageous from a processing efficiency 
point-of-vlew. 



According to yet another preferred embodiment of this aspect of 
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the invention, the metlioci involves an initial acquisition phase 
and a subsequent tracking phase. The acquisition phase esta- 
blishes a set of preliminary parameters which are required for 
initiating a decoding of signals that are received during the 

5 . tracking phase. The parameters may include: a modified code 
vector, a carrier frequency candidate vector an initial phase 
position, a code phase position and code index, which denotes a 
starting sample value for the modified code vector. A successful 
acquisition phase thus results in that at least one signal source 

10 specific code sequence is identified and that this signal is 
possible to track by the receiver (i.e. may be contihued to be 
received). , 

According to still another preferred embodiment of this aspect of 
the invention, the tracking In turn, involves the following steps. 

15 First, based the tracking characteristics, a prompt pointer is 
calculated for each modified code vector. The prompt pointer 
indicates a start position for the code sequence. The initial 
prompt pointer value is set equal to the code index. Then, at 
least one pair of early- and late pointers is assigned around 

20 each prompt pointer. The early pointer specifies a sample value 
which is positioned at least one element before the prompt 
pointer's position, and correspondingly, the late pointer specifies 
a sample value being positioned ait least one element after the 
prompt pointer's jaosition. These pointers are used to maintain 

25 the tracking of the received signal by a repeated repositioning of 
the pointers around a correlation maximum value, such that the 
. prompt pointer is positioned as close as possible to this value, 
the early pointer is positioned in time somewhat prior to this 
value and the late pointer is positioned in time somewhat after 

30 this value. An important advantage attained by this strategy Is 
that one and the same vector may be re-used in order to 
represent different delays. Thus, valuable memory space and 
processing capacity is sisived. 

According to still another preferred embodiment of this aspect of 
35 the invention, the tracking involves the following further steps. 
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First, a sequence of incoming level-discrete sample values is 
received. Then, data words are formed from the sample values, 
such that each data word contains a numbed of elements which 
is equal to the number of elements in each carrier frequency- 
5 phase candidate vector. After that, a relevant set of carrier 
frequency-phase candidate vectors is calculated for the data 
word, and a pre-generated in-phase representation respective a 
quadrature-phase representation of the vector is acquired for 
each carrier frequency-phase candidate vector in the relevant 

10 set. Subsequently, each data word is on one hand multiplied 
with a In-phase representation of the carrier frequency-phase 
candidate vector In the relevant set to produce a first inter- 
mediate-frequency-reduced information word, and on the other 
hand multiplied with a quadrature-phase representation of the 

16 carrier frequency-phase candidate vector in the relevant set to 
produce a second Intermediate-frequency-reduced information 
word. The intermediate-frequency-reduced information words 
thus constitute a respective component of the vectors mentioned 
above that have been derived from the incoming data words. 

20 This multiplied strategy is advantageous because it enables a 
very efficient usage of the capacity of a. general-purpose 
microprocessor. 

According to yet another preferred embodiment of this aspect of 
the invention, XOR- or SIMD-operations are used to perform the 
25 multiplication between the data word and the in-phase repre- 
sentation of the carrier frequency-phase candidate vector re- 
spective between the data word and the a quadrature-phase 
representation of the carrier frequency-phase candidate vector. 

According to another preferred embodiment of this aspect, of the 
30 invention, the tracking involves the further steps of: correlating 
the first intermediate-frequency-reduced information word with a 
modified code vector starting at a position Indicated by the 
prompt pointer to produce a first prompt-despread symbol string, 
correlating the first intermediate-frequency-reduced Informatiori 
36 word with a modified code vector starting at a position indicated 
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by the early pointer to produce a first early-despread symbol 
string, correlating the first intermediate-frequency-reduced Infor- 
mation word with a modified code vector starting at a position 
indicated by the late pointer to produce a first late-despread 

5 symbol string, correlating the second intermediate-frequency- 
reduced information word with a modified code vector starting at 
a position indicated by the prompt pointer to produce a second 
prompt-despread symbol string, correlating the second inter- 
mediate-frequency-reduced information word with a modified 

10 code vector starting at a position indicated by the early pointer 
to produce a second early-despread symbol string, and corre- 
lating the second intermediate-frequency-reduced information 
word with a modified code vector starting at a position indicated 
by the late pointer to produce a second late-despread symbol 

15 string. Preferably, a resulting data word is then derived for each 
set of the despread symbol strings. This may. for example be 
performed either by performing the relevant adding operation, or 
by looking up a pre-generated value In a table based on the bit 
pattern of the despread symbol strings depending on whether 

20 the resulting data words contain packed or un-packed 
information. Data of high quality may thus be obtained with a 
minimum of real time calculations, which is advantageous from a 
processing point-of-view. 

According to yet another preferred embodiment of this aspect of 
25 the invention, either XOR- or SIMD-operatlons are used also to 
perform the multiplication between the intermediate-frequency- 
reduced information words and the modified code vectors, since 
again, this is advantageous from a processing point-of-view. 

According to another preferred embodiment of this aspect of the 
30 invention, certain pieces of information are propagated in 
connection with completing the processing of a current data 
word and initiating the processing of a subsequent data word. 
Preferably, this propagated information includes: a pointer which 
indicates a first sample value of the following data word, a group 
35 of parameters that describe the relevant set of carrier 
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frequency-phase candidate vectors, the relevant set of code 
vectors, and the prompt-, early-, and late pointers. Hence, based 
on the propagated information, the subsequent data word may 
be processed immediately. 

5 According to a further aspect of the invention the object is 
achieved by a computer program directly loadable into the 
internal memory of a computer, comprising software for perfor- 
ming the above proposed method wheri said program is run on a 
computer. 

10 According to another aspect of the invention the object Is 
achieved by a computer readable medium, having a program 
recorded thereon, wliere the program is to malce a computer 
perform the above proposed method. 

According to another aspect of the invention the object is 
15 achieved by the signal receiver as initially described, which is 
characterized in that the proposed coniputer program is loaded 
Into the memory means, such that the receiver will operate 
according to the above-described method. 

Consequently, the invention offers an excellent instrument for 
20 receiving and decoding any type of spread spectrum signals in a 
general software radio receiver. For instance, one and the same 
receiver may be used to tracic navigation satellite signals of 
different formats. These signals may even be received in 
parallel, such that a first signal is received on a first format 
25 while one or more other signals are received on a second or 
third format. 

Moreover, the proposed solution permits one and the same 
receiver to be utilized for completely different purposes, such as 
communication in a cellular system. 

30 BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention is now to be explained more closely by 
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means of preferred embodinnents. which are disclosed as 
examples, and with reference to the attached drawings. 

Figure 1 shows a spectrum graph over an exemplary spread 
spectrum signal from which Information may be ex- 
5 tracted according to the invention, 

Figure 2 shows a spectrum graph over a frequency down- 
converted signal corresponding to the spectrum 
shown in figure 1 , 

Figures is an enlarged version of figure 2 in which a 
10 Doppler shifting of the spectrum is illustrated, 

Figure 4 illustrates how a data signal is modulated onto a 
signal source specific code sequence according to 
an embodiment of the invention, 

Figure 5 depicts a three-dimensional graph over carrier 
15 frequency-phase candidate vectors to be used as 

initial correlators according to an embodiment of 
the invention. 

Figure 6 . shows an aiteirnatiye representation the carrier 
frequency-phase candidate vectors in figure 5, 

20 Figures 7r8 show graphs which illustrate how a code index and 

a code phase are defined for a received signal 
according to an embodiment of the invention. 
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Figures 9a,b illustrate how a modified code vector is generated 
based on an original code vector according to an 
embodiment of the Invention, 



Figure 10 



30 Figure 11 



depicts a three-dimerisional graph over code rate- 
phase candidate vectors to be used as subsequent 
correlators according to an embodiment of the 
invention. 

illustrates how a set of code sequence start 
pointers is assigned to a modified code vector 
according to an embodiment of the invention. 
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Figure 12 
Figure 13 

5 

Figure 14 

10 Figure 15 
Figure 16 



illustrates how data words are formed from a 
sequence of incoming level-discrete sample values 
according to an embodiment of the invention, 

illustrates how the data words of figure 12 are 
correlated with various pre-generated vectors 
according to an embodiment of the invention, 

illustrates how multiple pairs of pointers may be 
assigned according to an embodiment of the 
invention, . 

shows a proposed signal receiver, and 

illustrates, by means of a flow diagram, the general 
method of processing spread spectrum signals 
according to the invention. 
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DESCRIPTION OF PREFERRED EMBODIMENTS OF THE 
15 INVENTION 

In order to realize a truly software based signal processing 
solution in real-time at high sampling rates, such as those 
required in a GNSS, the architecture of the microprocessor 
system must be utilized very efficiently. The present invention 

20 uses an approach, u^ich involves a form of parallel processing 
based on pre-generated code vector representations. The 
invention thereby strikes a balance between processing speed 
and memory usage, which compensates for the capacity deficit 
of a general purpose microprocessor in relation to an ASIC. The 

25 details of this solution will become apparent from the below 
disclosure. 

Figure 1 shows a spectrum graph over an exemplary spread 
spectrum signal from which information may be extracted 
according to the invention. For example, the signal may be a 
30 GPS C/A-signai in the LI -band with a carrier frequency fnp of 
1.57542 GHz. This means that the signal's spectrum is 
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symmetrical around the carrier frequency fHF- Moreover, the 
signal has at least 98% of its energy distributed within a 2 MHz 
wide frequency band BWhf around the carrier frequency fHF- 
Therefore, it is generally sufficient If the spread spectrum signal 
5 is downconverted to an intermediate frequency fur at around 1 
MHz In order to enable further signal processing at a lower 
frequency.. 

Figure 2 shows a spectrum graph over such a frequency down- 
converted signal, whose spectrum Is symmetrical around an 

10 intermediate frequency fip. which Is considerably lower than the 
carrier frequency fHF, for example fiF = 1 ,25 MHz. Due to move- 
ments of the transmitter (i.e. typically a satellite in orbit) and the 
receiver (i.e. typically a mobile GNSS-receiver) a Doppler shift 
may occur in the received signal, which is proportional to the 

15 relative velocity between the transmitter and the receiver. A 
maximum Doppler shift fo Is defined which In the GNSS case 
normally is approximately 10 kHz. Thus, the entire spectrum 
may be shifted within an interval 2fD, such that the center 
. frequency falls any where between fiF+fo (corresponding to a 

20 maximum positive Doppler shift) and fip-fo (corresponding to a 
maximum negative Doppler shift). For illustrative purposes the 
Interval 2fo. has been severely exaggerated in the figure 2. The 
. figure also indicates a phase variation interval <pd (along an 
independent axis <p) around the Intermediate frequency fiF. The 

25 phase variation interval 90 demonstrates the fact that the phase 
position of the received signal is initially unknown to the 
receiver, and thus constitutes a calibration parameter. The 
effects of this parameter will be further elucidated below with 
reference to figure 5. 

30 An enlarged version of the spectrum In figure 2 is shown in 
figure 3. The intermediate spectrum may, due to possible 
Doppler effects, be shifted such that its center frequency lies in 
the interval 2fo from f|F.min=fiF-fo to fiF-max^fiF+fo- The potential 
maximum shifted spectra are here indicated with a respective 

35 dashed line. The diagram in figure 3 also shows a Nyquist 




16 



frequency rs/2 corresponding to a minimal sampling rate rs 
which, when used for sampling the intermediate frequency 
signal, will result in a non-aliased discrete spectrum. Provided 
that fo » 10 Miz, fiF s 1,25 MHz. BWhp » 2 MHz and the inter- 
5 mediate frequency signal is lowpass filtered through an analog 
filter having a 2,5 MHz wide passband, the minimal sampling 
rate rs becomes 5 MHz (i.e. two times the highest frequency 
component after the analog filter). 

Figure 4 illustrates schematically how a data signal D is 

10 modulated onto a signal source specific code sequence CS on 
tiie transmitter side according to an embodiment of the inven- 
tion. According to a preferred embodiment of the invention, the 
signal source specific code sequence CS constitutes so-called 
pseudo random noise. The data signal D here contains a data 

15 symbol sequence [-i-l, -1, -1, +1, -1], and has a relatively low 
symbol rate, say 50 Hz. The signal source specific code 
sequence CS, however, has a relatively high symbol rate (or 
more correctly chipping rate). For instance, a signal source 
specific code sequence CS in the form of a GPS C/A code may 

20 have a chipping rate of 1,023 MHz and contain 1023 chips per 
period. Each chip, ch, is either +1 or -1. Hence, the C/A code 
repeats itself one per ms. The data signal D is modulated onto 
(or spread by) the signal source specific code sequence CS by 
multiplying each data symbol with the code sequence. CS. The 

25 data symbol -<-1 thereby results in an unaltered code sequence 
CS, whereas the data symbol -1 results in an inverted code 
sequence CS. 

If a data signal D having a rate of 50 Hz is spread by means of a 
signal source specific code sequence CS having a chipping rate 
30 of 1,023 MHz, this results in 20 entire code sequences CS per 
data symbol. Namely, the period time for one data symbol is 20 
ms, whereas the period time for the code sequence is only 1 ms. 

Figure 5 depicts a three-dimensional graph over different carrier 
frequency-phase candidate vectors, which may be used as initial 
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correlators according to an embodiment . of the invention. As 
mentioned above, the intermediate frequency spectrum may be 
Doppler shifted. Therefore, the particular frequency of the 
received signal is initiaily unl(nown within a Doppler frequency 

5 interval from fiF-mm to f|F-max- (Moreover, the phase position of the 
received signal is beforehand unl^nown to the receiver. In order 
to enable an estimation of the actual intermediate carrier fre- 
quency, the Doppler frequency interval fiF^nin ' fiF-max ie divided 
Into an integer number of equidistant frequency steps. Of course, 

10 the larger the number of steps is, the more accurately the inter- 
mediate carrier frequency may be estimated. In this example, a 
step size Af of 20 Hz is chosen. In some applications, however, a 
much larger step size Af may provide a sufficient frequency 
accuracy, whereas in other applications (where relatively high 

15 accuracy is needed) a considerably smaller step size Af may be 
required. In any case, a Doppler frequency interval fiFLmin - ^iF^max of 
20 kHz and a step size Af of 20 Hz requires 1000 steps, where a 
first step includes the frequencies from ftF-mm to fiF.min+Af, a second 
step includes the frequencies from fiF.min+Af to fiF.min+2Af, and so 

20 on up to a last step including Uie frequencies from fiF.inin-''999Af to 

f|F-max' 

A frequency candidate vector fiF.c is defined for each frequency 
step. In order to estimate an initial phase position of the 
intermediate carrier frequency, an integer number of initial phase 

25 positions %, <i>7 are determined. The illustrated example 
shows eight such positions <Pc. This gives a phase resolution of 
n/4. For each frequency candidate vector fiF^, eight different 
phase positions %, 9/ are thus defined which each repre- 
sents a particular phase shift from 0 to JnlA. A combination of a 

30 certain carrier frequency candidate vector fiF-c and a particular 
initial phase position 9c is referred to as a carrier frequency- 
phase candidate vector. Figure 5 shows all such vectors (i.e. 
eight) for the first step, where ftp^ ~ fie^m, end the second step, 
where fiF^ « fiFnnm+Af. 

35 Additionally, figure 5 includes an axis Lc along which the 
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number of sample values (or length) of the frequency candidate 
vector fiF^ is indicated. The Lc number is preferably a power of 2 
because this generally renders the implementation more simple. 
The number of sample values in each frequency candidate vector 

5 fiF.c is set in consideration of the available memory space, the 
desired frequency resolution and the desired phase resolution. If, 
in the example above, a length of 512 samples is chosen and each 
sample requires one byte, the total (uncompressed) memory 
demand would be 3,9 MB (i.e. 1000 frequencies x 8 phase 

10 positions x 512 B » 4098000 B » 3,9 MB). 

Figure 6 shows an alternative representation the carrier 
frequency-phase candidate vectors Vf,(fiF-c. *c) in figure 5. 
. Here, the vectors Vf,(f|F.c, <Pc) range from the carrier frequency 
fiF-c=fiF.min (say 1.24 MHz) to fiF-c=fiF.max (say 1.26 MHz), from 
15 the phase shift position <Pc«<Po (say 0) to 9c=<p7 (say 7«/4). all 
have a length Lc of s„ samples (say 512) and are represented in 
a cube. Hence, a particular frequency-phase candidate vector 
Vfq,(f|F.c* ^c) is given by a cuboid segment in the cube which is 
parallel with the Lc~axis. and ranges from Lc=si to Lc~Sn. 

20 . After having generated all the carrier frequency-phase candidate 
vectors Vfp(f|F.Cf *Pc). the vectors are stored in a digital memory 
with a comparatively short access time, such that they may be 
readily accessed during a later acquisition and/or tracking 
phase. Preferably, the carrier frequency-phase candidate 

25 vectors Vf9(f|F.ci ^c) are stored according to a data format which 
is adapted to performing correlation operations between a vector 
that has been derived from incoming data words (representing 
one or more received signals) and a segment of a carrier 
frequency-phase candidate vector Vf^(fiF-c. ^c)- The vectors may 

30 either include packed sampies or they may contain samples that 
are unpacked to a smallest manageable unit in the micro- 
processor, such as a byte. Further details regarding the data 
format Issues will be discussed below. 

Figure 7 shows a graph over an exemplary code sequence CS, 
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which has been received, downcon verted and sampled 
according to an embodiment of the invention. After subsequent 
removal of the intermediate frequency cornponent, the code 
sequence CS either attains the signal value -1-I or -1, and the 

5 symbols are altered at a chipping rate of for example 1,023 
MHz. At the receiver side, the code sequence CS is sampled 
with a sampling rate rs of for example 5 MHz. The sampling in- 
stances are indicated by means of dots on the code sequence 
signal CS along an axis S|. As is apparent from the figure there 

10 are approxirhately four sample values per chip period. Hence, 
the requirements of the Nyquist theorem are clearly fulfilled. 
However, the samples are not capable of determining the 
transition Instances with a sufficient accuracy, i.e. where the 
code sequence signal CS changes from representing one chip 

15 symbol to representing another one chip' symbol. Therefore, a 
so-called code phase must be determined. Moreover, a parti- 
cular sample value at which the code sequence starts must also 
be established. 

Figure 8 shows a graph, which illustrates the sample value at 
20 which the code sequence CS starts (the so-called code index 
CI) and how a code phase, Cph, is defined in relation to this 
sample value. An initial segment of the code sequence signal 
CS in figure 7 is therefore represented in figure 8. 

A first sample instance is here presumed to occur shortly prior 
25 to the actual start of the code sequence CS signal period. 
Hence, a corresponding sample value Si belongs to the end of 
the preceding period. This is illustrated by means of the line 
representing the code sequence signal CS being dashed. A 
subsequent sample instance, however, overlaps the code 
30 sequence CS period. The corresponding sample value sa is 
defined as the code index CI. The actual beginning of a code 
sequence signal is determined in a preceding acquisition phase 
and typically involves correlation between the received signal 
and a local copy of the code sequence CS. 




• « • ■ • « • 
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The distance between the start of the code sequence signal CS 
and the code index Ct. in turn, is defined as the code phase 
Cph. land thus represents a measure of the translation (or skew) 
between the sample instances Si and the chip transitions. The 

5 sampling interval between two consecutive sample instances Si 
is divided into an integer number of possible Initial code phase 
positions in order to estimate the code phase Cph. In this 
example, ten 0.1 -steps ranging from 0.0 to 0.9 are defined, 
where a first step 0.0 would indicate that code sequence signal 

10 CS starts at the code index CI, and a last step 0.9 would 
indicate that code sequence signal CS starts aiiriost at the 
preceding sample value Si. The figure 8 shows a code phase 
Cph of 0.4, i.e. a situation where the code sequence signal CS 
starts approximately in the middle between two consecutive 

15 sample instances S{. 

In analogy with the carrier frequency-phase candidate vectors 
■ Vf9(fiF.c. <Pc), code rate-phase candidate vectors may now be 
defined for each signal source specific code sequence CS and 
combination of code rate CR. code Index CI and code phase 
20 Cph. Nevertheless, according to a preferred embodiment of the 
invention, a modified code vector is first produced based on 
each original code vector. 

Figures 9a and 9b illustrate how this may be performed. A code 
vector CV in figure 9a contains a number of sample values, say 

25 5000. Preferably, the number of sample values is selected in 
relation to sampling rate such that a full code sequence period 
Is recorded, in the above example, with a sampling rate rs of 5 
MHz, this means that the 5000 samples cover exactly one full 
code sequence CS period (which has a duration of 1 ms and 

30 contains 1023 chips, since the chipping rate is 1,023 MHz). 

A modified code vector CVm is now generated on basis of the 
original code vector CV by copying a particular number of 
elements Ee (e.g. two) from the end of the original code vector 
CV to the beginning of the modified code vector CV^. 
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Correspondingly, a particular number of elements Eb (e.g. two) 
from the beginning of the original code vector CV are copied to 
the end of the code vector . CVm. i.e. after the elements Ee. 
Figure 9b shows the resulting modified code vector CVm. Since. 

5 in this example, each chip symbol is sampled at approximately 
four instances, two sample values are equivalent to about one 
half chip period. Thereby, the above copying of the elements E© 
and Eb allows a one-haif-chip shifting of the code vector CV 
during a correlation operation between the code vector CV and 

10 the received data. This will be discussed further below with 
reference to figure 11. 

Figure 10 depicts a three-dimensional graph over different code 
rate-phase candidate vectors CVm(CRc.c. Cph) which represent 
a signal source specific code sequence CS(i). Analogous to the 
15 intermediate carrier frequency, the nominal code rate CRc for 
the code rate candidate vector CRc-c may also vary due to 
Doppler effects from a lowest value CRmm to a highest value 
CRmax- However, since the code rate (or chipping rate) is com- 
paratively high, the rate variation here becomes rather small. 

20 For example, a Doppler shift of 10 icHz at the intermediate 
carrier frequency is equivalent to a shift in a nominal code rate 
CRc at 1,023 MHz of merely 6.5 Hz. Namely, the Doppler shift In 
respect of the code rate is generally equal to the ratio between 
the carrier frequency and the code rate (1575.42 MHz / 1,023 

25 MHz » 1540; a Doppler shift of 10 kHz at the intermediate 
carrier frequency thus corresponds to 10000/1540 « 6.5 Hz). 
The code rate candidate vector CRc-c may attain a value any- 
where in a Doppler rate Interval between CRmm - 1022993,5 Hz 
and CRmax ~ 102306,5 Hz. The Doppler rate interval CRmm - 

30 CRmax is divided into an integer number of equidistant code rate 
steps ACR, say 13. In this example, each interval CRmm. 
CRmm'i'ACR. .... CRmax thus spans 1 Hz. According to a preferred 
embodiment of the invention, the code rate step size is adaptive, 
e.g. determined on basis of the sampling frequency and the 

35 desired resolution. 
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In figure 10, the three-dimensional graph (or cube) contains a set 
of code rate-phase candidate vectors CVm(CRc.c. Cph). where the 
code rate ranges from CRc^sCRmin to CRc^'^CRmax. the code 
phase ranges from Cph«0.0 to CphsO.9. and each vector has a 
5 length Lcvm of Sm samples (for example 5004). Hence, a particular 
code rate-phase candidate vector CVm(CRc.c. Cph) is given by a 
cuboid segment CVm(CRn,in+2ACR, 0.4) which is parallel with the 
Lcvm -axis and ranges from Lcvm-^i to Lcvm =^3^. 

After having generated all the code rate-phase candidate 
10 vectors CVm(CRc.c. Cph), these vectors are stored in a digital 
memory with a comparatively short access time, such that they 
may be readily accessed during the acquisition and/or tracking 
phase. Preferably, the code rate-phase candidate vectors 
CVm(CRc^, Cph) are stored according to a data format which is 
15 adapted to correlation operations to be performed between a 
vector that has been derived from incoming data words and a 
segment of a the code rate-phase candidate vectors CV„(CRc.c, 
Cph). For example, the code rate-phase candidate vectors 
CV,n(CRc.c. Cph) may include packed samples or they may 
20 contain samples that are unpacked to a smallest manageable 
unit in the microprocessor, such as a byte. Further details 
regarding these data format issues will be discussed below with 
reference to figure 13. 

During tracking of a spread spectrum signal, i.e. a continued 
25 reception of the signal and demodulation of the data contained 
therein, it is necessary to update relevant tracking parameters 
for the signal. Normally, the tracking phase Is preceded by a so- 
called acquisition phase during which a set of preliminary 
parameters are established that are required for initiating the 
30 signal decoding. A successful acquisition phase thus results in 
that at least one signal source specific code sequence is 
identified and a data signal transmitted by means of this 
sequence Is possible to demodulate. The tracking characteris- 
tics that may be associated with a signal source specific code 
35 sequence include: a modified code vector CVm (which defines a 
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carrier frequency candidate vector fip^; and an initial phase 
position <l>c). a code phase position, Cpli, and a code index, CI. 
In order to maintain an adequate timing of the correlating 
operations in the receiver, a set of pointers are required which 
5 indicate the start position of the code sequence in relation to the 
modified code vector and are updated repeatedly, preferably 
between each correlation. 

Figure 1 1 illustrates one such set of pointers that are assigned 
to a modified code vector CVm. A prompt pointer Pp indicates a 

10 current estimate of the code sequence start position. This 
pointer Pp is initially set equal to the code index CI. Adequate 
Pp-pointer positions for any subsequent segments of the 
received data are obtained by retrieving appropriate modified 
code vectors C\/„ from the digital memory where such pre- 

15 generated are stored. Between each segment, the code rate 
candidate vector CRc-c as well as the initial code phase Cph is 
updated. Additionally, at least one pair of early- and late 
pointers Pe and P|. respectively is assigned on each side of the 
prompt pointer Pp. where the early pointer Pe specifies a sample. , 

20 value being. positioned at least one element before the prompt 
pointer's Pp position, and the late pointer Ft specifies a sample 
value being positioned at least one element after the prompt 
pointer's Pp position. Each correlation is then performed with the 
presumption that the code sequence starts at the Pp-position, at 

25 the PE-position as well as at the Pt-position. If the correlation 
over the Pp-position results in a higher correlation value than 
over any of the other positions, and at the same time the Pe- and 
Pt-pointers are balanced (i.e. the correlation over these 
positions result in equal values), this is interpreted as an optimal 

30 setting of the tracking parameters. Should, however, this not be 
true the sampled data must be repositioned in relation to the Pp-. 
Pe- and Pt-pointers, such that the correlation peak again 
coincides with the Pp-pointer's position. An un-callbrated pointer 
positioning may thus be detected via a difference between the 

35 correlation value over the Pe-pointer position and the correlation 

i 
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value over the PL-pointer position, while the correlation over the 
Pp-position still results in the highest correlation value. In this 
case, the potentially highest correlation value (corresponding to 
the optimal Pp-positioning) lies somewhere between the current 
5 Pp-position and the current Pe-position or between the current 
Pp-position and the current Pc-posltion depending on which of 
the Pe- and Pt-pointer positions that results in the highest 
correlation value. A pointer updating normally takes place 
between each segment. A second or third pair of Pe- and Pl- 

10 pointers outside the first pair of pointers further enhances the 
possibilities of accomplishing an accurate adjustment of the 
tracicing. Depending on the value of the received symbol, an 
optimal setting Pp-pointer may, in fact, also be equivalent to the 
correlation over the Pp-position resulting in a lower correlation 

15 value than over any of the other pointer positions. Namely, this 
is true if a negative valued data bit is received. However 
otherwise, the same principles apply. 

In any case, the Pp-, Pe- and pL-pointers may only be positioned 
within certain intervals in the modified code vector CVm. The Pp- 

20 pointer may. be set anywhere within a first interval 111, the Pg- 
. pointer within a second interval 112 and the Pt-pointer within a 
third interval 113. Ail the pointers are always set together, such 
that their relative distances remain unchanged. A limit value PL-min 
for the Pt-pointer typically coincides with the first element of the 

25 modified code vector CVm. Namely, outside this range a 
meaningful correlation cannot be performed. As a consequence of 
this and the fact that the pointers are set jointly Pp-, Pe- and Pl, a 
corresponding limit value Pe^jn for the Pe-pointer is also defined. 

Figure 12 illustrates how data words d(1). d(2), .... d(N) are 
30 formed from a sequence 1210 of incoming level-discrete sample 
values according to an embodiment of the invention. The data 
words d(1). d(N) are formed such that each data word d(k) 
contains a number of elements which is equal to the number of 
elements Sn in each carrier frequency-phase candidate vector. For 
35 instance, the length of these vectors may be 512 sample values 
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(as In the example described above with reference to the figures 
5 and 6). Given that each sequence 1210 contains 5000 sample 
values, the number N of data words becomes 10 (5000/51 2»9,76 
means that the last quarter of the tenth data word d(10) Is empty). 

5 According to the invention, the elements in a current data word 
d(k} are processed jointly (I.e. in parallel). Then, the processing of 
a subsequent data word d(k-t-l) Is Initiated, and any processing 
parameters obtained during the processing of the current data 
word (i{k) that are required when processing the subsequent data 

10 word d(i(-t-i) are propagated as input data to the processing of the 
latter. These processing parameters preferably include: a pointer 
Pd which Indicates a first sample value of the subsequent data 
word dCk+l), a group of parameters which describe the relevant 
set of carrier frequency-phase candidate vectors Vf,,(fiF.c, ^c) (i-e* 

15 a carrier frequency candidate vector fiF.c representing an in-phase 
version and a carrier frequency candidate vector fip.c representing 
a quadrature-phase version), the relevant set of code vectors 
CVm, and prompt-, early-, and late pointers Pp, Pe. and Pl 
respectively. 

20 Figure 13 illustrates how the data words of figure 12 are corre- 
lated with various pre-generated vectors according to an embo- 
diment of the invention. First, a relevant set of carrier frequency- 
phase candidate vectors Vf<p(f|p^, <Pc) is calculated for the data 
word d(k). This set is normally obtained In a preceding 

25 acquisition phase (or, if the data word d(k) is not the first word 
to be processed, during the processing of a previous word). 

Thereafter, for each carrier frequency-phase candidate vector 
Vfq.(fiF.c. ^c) in the relevant set, a pre-generated in-phase repre- 
sentation fiFi and a quadrature-phase representation fipo of the 

30 vector (Vf^(fiF.c. ^c)) respectively is acquired. Both these repre- 
sentations are found In the digital memory where the pre-gene- 
rated carrier frequency-phase candidate vectors Vf,(fiF.c, ^c) are 
stored because they are merely one quarter period delays of 
one and the same vector, and the memory holds initial phase 

35 position representations for an entire period, see figure 6. 
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Then, in order to eliminate the influence of the carrier frequency 
component, the data word d{k) is multiplied with the in-phase 
representation fiFi of the carrier frequency-phase candidate 
vector Vf,(fiF.c> ^^c) in the relevant set. As a result thereof, a first 
intermediate-frequency-reduced information word SiF.i(k) is pro- 
duced. Correspondingly, the data word d<k) is also multiplied 
with the quadrature-phase representation fipo of the carrier 
frequency-phase candidate vector Vf^(fiF-c, ^c) «n the relevant 
set, and a second intermediate-frequency-reduced information 
word S|f.q(I<) is produced. 

Preferably, the multiplication between the data word d(l() and 
the in-phase representation fiFi respective the quadrature-phase 
representation fipQ of the carrier frequency-phase candidate 
vector Vfp(fiF.c. ?c) are performed by means of XOR- or SII^O- 
operations (I.e. 1-bit and multiple-bits multiplications respec- 
tively). This is namely possible If the data word d(l() and the 
vector Vf,(fiF^, <Pc) have compatible data formats. 

Subsequently, the first intermediate-frequency-reduced informa- 
tion word SiF.i(k) is multiplied with a modified code vector CVm.p, 
which has been retrieved from the digital memory where these 
pre-generated vectors are stored, and starts at a position 
indicated by the prompt pointer Pp. A first prompt-despread 
symbol string Aip(k) is produced as a result of this operation. 
The first intermediate-frequency-reduced information word 
SjF.i(k) is also multiplied with a modified code vector CVm.E(k). 
which starts at a jposition indicated by the early pointer Pg.. A 
first eariy-despread symbol string AiE(k) is produced as a result 
of this operation. The first intermediate-frequency-reduced 
information word SiF-i(k) is likewise multiplied with a modified 
code vector CVm.L(k), which starts at a position indicated by the 
late pointer Pl, and a first late-despread symbol string AiL(k) is 
produced. Moreover, the second intermediate-frequency- 
reduced information word S,F.Q(k) is multiplied with a modified 
code vector CVm,p(k), which starts at a position indicated by the 
prompt pointer Pp, and a second prompt-despread symbol string 
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AQp(k) is produced. Correspondingly, the second intermediate- 
frequency-reduced information word Sif.q(Ic) is multiplied with a 
modified code vector CVm.E(k), which starts at a position 
indicated by the early pointer Pe. and a second eariy-despread 
5 symbol string AqeCIc) is produced. Liicewise, the second 
intermediate-frequency-reduced information word Sif.q(<c) is 
multiplied with a modified code vector CVm.L(l<)> which starts at a 
position indicated by the late pointer Pl, and a. second late- 
despread symbol string AQL(k) is produced. 

10 After that, a resulting data word DR.iE(k), OR.ip(k), DR.iL(k), 
DR.Qe(k), DR.Qp(k) and DR.QL(k) may be derived from the each of 
the despread symbol strings Aip(k). AiE(k), AiL(k), AQp(k), AQE(k) 
and AQL(k) respectively by adding the elements in the respective 
string together. If the strings represent un-packed data, the 

15 processor may simply perform the relevant adding operation to 
obtain the resulting data word DR(k). If however, the strings 
represent packed data, the resulting data word DR(k) is, 
according to a preferred embodiment of the invention, derived by 
looking up a pre-generated value in a table 1310 based on the 

20 bit pattern given by the respective despread symbol string 
A,p(k). AiE(k), A,L(k), AQp(k), AQE(k) and AQ^k). 

According to a preferred embodiment of the invention, in addition 
to the individual data words DR.iE(k), DR.|p(k), DR.iL(k), DR.QE(k), 
DR>Qp(k) and DR.QL(k), corresponding accumulated data words are 
25 produced, such that after having generated a data word OR(k), the 
sum of all data words Dr(1) to OR(k) is also obtained. I4ence, 
when the data word d(N) has been processed, the resulting sum 
Dr(1 ) to Dr(N) is also at hand. 

Moreover, the resulting data words represent payload information. 
30 For Instance, based on the set of resulting data words DR.iE(k), 
DR-ip(k), DR^t(k). DR^E(k), DR.Qp(k) and DR.^L(k) a demodulated 
piece of information of a transmitted data symbol sequence D as 
shown In figure 4 may be derived. In a steady state operation, 
basically only one of the resulting data words, say OR.rp(k), 
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actually carries the information. 

According to a preferred embodiment of the invention, the multi- 
plications involved in the multiplications between the interme- 
diate-frequency-reduced information words S|(:.|(l() and S|f:.Q(k) 
5 respectively and the modified code vectors CVm.p(k), CVm.E(k}; 
CVro.i.(k) are all performed by means of XOR- or SiMD-ope- 
rations. This Is possible if the intermediate frequency infor- 
mation words SiF-i(l«); SiF.Q(k) and the modified code vectors 
CVm.E(ic). CVm.p(k}; CVm.L(K) have compatible data formats. 

10 It was mentioned above with reference to figure 11 that more 
than one pair of early- and late pointers may be assigned to the 
prompt pointer in order to sidcompiish an improved tracking. 
Figure 14 illustrates how two pairs of such pointers Pei; Pli and 
Peal Pl2 respectively are assigned around a prompt pointer Pp 

15 according to an embodiment of the invention. The horizontal 
axis here represents a code shift CS and the vertical axis shows 
a normalised correlation factor. An optimal positioning of the 
prompt pointer Pp is thus equivalent to a normalised correlation 
value of 1 . The more the sampled data is shifted in either 

20 . direction from this position, the lower the correlation value 
becomes down to substantially zero for a shift of one chip or 
more. The algorithm is perfectly balanced when (as illustrated in 
the figure) both the pointers Pei and Pli of a first pair produce a 
first shifted correlation result, both the pointers Pe2 and P1.2 of a 

25 second pair produce a second shifted correlation result, the first 
shifted correlation result is larger than the second shifted corre- 
lation result and both the first and the second shifted correlation 
results are smaller than the correlation result obtained at the 
prompt position Pp. Said correlation results are preferably 

30 generated according to a process corresponding to what has 
been described above with reference to figure 13. Based on the 
correlation results for the pointer positions Pp. Pei, Pli> Pe2 and 
Pl2 a (if necessary new) set of modified code vectors CVm.E(k), 
CVm.p(k): CVn,.L(k) is selected, which is deemed to be optimal. 
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Figure 15 shows a signal receiver 1500 for receiving navigation 
data signals transmitted in a navigation satellite system accor- 
ding to an embodiment of the invention. The receiver 1500 
includes radio front end unit 1510, an interface unit 1520 and a 
5 digital processor unit 1530. 

The radio front end unit 1510 is adapted to receive a continuous 
radio signal Shf. and in resporise thereto produce a corre- 
sponding electrical signal S|f which has a comparatively high 
frequency. The interface unit 1520 is adapted to receive the 

10 electrical signal S|(r, and in response thereto, produce a 
sequence of sample values that represents the same information 
as the electrical signal Sif and is divided into data words d{k). 
. The digital processor unit 1530 is adapted to receive the data 
words d(lc), and in response thereto, demodulate a data signal. 

15 The digital processor unit 1530, in turn, includes a memory 
means 1535, which is loaded with a computer program that is 
ciapable of controlling the steps of the proposed procedure when 
the program is run in the processor unit 1530. 

In order to sum up, the general method of processing spread 
20 spectrum signals according to the invention will now be 
described with reference to a flow diagram in figure 16. 

A preparation step 1600 pre-generates code vectors which each 
represents a signal source specific code sequence that is 
intended to be received (or at least possible to receive) and 
25 demodulated in the receiver. The step 1600 is performed before 
the signal reception is initiated. 

A step 1610 then receives a continuous signal of a compara- 
tively high frequency. A following step 1620 samples the 
continuous signal at a basic sampling rate, whereby a resulting 
30 sequence of time discrete signal samples is produced. Each 
sample is also quantised (either with a relatively low-resolution, 
such as with 1 bit per sample value, or with a relatively high 
resolution depending on the application and the number of data 
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bits per samples delivered from the radio front end unit 1510), 
such that a corresponding level-discrete sample value is 
obtained. Subsequently, a step 1630 forms data words from the 
sample values, where each data word includes one or more 
5 consecutive sample values. After that, a step 1640 carries out 
correlation operations between the information in the data words 
and at least one of the pre-generated representations of a signal 
source specific code sequence. The correlation step involves 
correlating at least each vector in a sub-group of the code 

10 vectors with at least brie vector that lias, been derived from a 
data word. A following step 1650 produces data as a result of 
the correlation performed in step 1640. Then, a step 1660 
investigates whether the sampled sequence has ended, i.e. 
whether there are any more data words to process. If more 

15 sampled data to process is found, the procedure returns to the 
step 1640. Otherwise, the procedure loops back to the step 
1610 again for continued reception of the incoming signal. 
Naturally, such signal is preferably also received during 
execution of the steps 1620 to 1660. The sequential procedure 

20 illustrated in figure 16 is merely applicable to a specific received 
signal segment. Preferably, all the steps are actually 1610 to 
1660 are actually performed in parallel. 

The process steps, as well as any sub-sequence of steps, 
described witii reference to the figure 16 above may be 

25 controlled by means of a programmed computer apparatus, such 
as a microprocessor located in a GNSS-receiver. Moreover, 
although the embodiments of the invention described above with 
reference to the drawings comprise computer apparatus and 
processes performed in computer apparatus, the invention thus 

30 also extends to computer programs, particularly computer 
programs on or In a carrier, adapted for putting the invention 
into practice. The program may be in the form of source code, 
object code, a code intermediate source and object code such 
as in partially compiled form, or in any other form suitable for 

35 use in the implementation of the process according to the 
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invention. The carrier may be any entity or device capable of 
carrying tlie program. For example, the carrier may comprise a 
storage medium, such as a ROM (Read Only Memory), for 
example a CD (Compact Disc) or a semiconductor ROM, or a 

5 magnetic recording medium, for example a floppy disc or hard 
disc. Further, the carrier may be a transmissible carrier such as 
an electrical or optical signal which may be conveyed via 
electrical or optical cable or by radio or by other means. When 
the program is embodied in a signal which may be conveyed 

10 directly by a cable or other device or means, the carrier may be 
constituted by such cable or. device or means. Alternatively, the 
. carrier may be an integrated circuit in which the program Is 
embedded, the integrated circuit being adapted for performing, 
or for use in the performance of, the relevant processes. 

15 The term "comprises/comprising" when used in this specification 
is taken to specify the presence of stated features, integers, 
steps or components. However, the term does not preclude the 
presence or addition of one or more additional features, 
integers, steps or components or groups thereof. 

20 The invention is not restricted to the descrifa|ed embodiments In 
the figures, but may be varied freely within the scope of the 
claims. 
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Claims 

1. A method of processing spread spectrum signals com- 
prising: 

receiving a continuous signal (S|f) of a comparatively higli 
5 frequency. 

sampling the continuous signal (Sif) at a basic sampling rate 
(rs), whereby a resulting sequence of time discrete signal 
samples (S[s;]) is produced, 

quantising each signal sample into a corresponding levet- 
10 discrete sample value, 

forming of a plurality of data words (d(1), .... d(N)) which 
each includes one or more consecutive sample values (Si, .... Sn). 
and 

correlation between information in the data words (d(k)) and 
15 at least one representation (CS(i)) of a signal source specific 
code sequence (CS), characterized by 

the method comprising a preparation for the correlation step, 
wherein, before receiving the continuous signal <Sif). a multitude 
of code vectors (CV; CVm) are pre-generated. each code vector 

20 representing, a particular code sequence (CS(i)) of the at least 
one signal source specific code sequence (CS), and 

the correlation step involving multiplying at least each vector 
in a sub-group (CVm.E, CVm.p; CVm.|.) of the code vectors (CV; 
CVm) with at least one vector (SiF.i(k): S|F^(k)) derived from the 

25 data word (d(ic». 

2. A method according to claim 1, characterized by each 
code vector (CV; CVm) representing a particular signal source 
specific code sequence (CS(i)) being sampled at the basic sam- 
pling rate (r,) and quantised with the quantising process being 

30 used to produce the level-discrete sample values (si. .... s„). 

3. A method according to any one of the preceding claims, 
characterized by the at least one signal source specific code 
sequence (CS) being pseudo random noise. 
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4. A method according to any one of the preceding claims, 
characterized by the receiving step involving down conversion 
of an incoming high-frequency signal (Shf) to an Intermediate 
frequency signal (Sif), the high-frequency signal (Shf) having a 

5 spectrum which is symmetric around a first frequency (fw) and 
the intermediate frequency signal (Sif) having a spectrum which 
is symmetric around a second frequency (fip) being considerably 
lower than the first frequency (fHF)- 

5. A method according to claim 4. characterized by the steps 
10 of: determining a maximum frequency variation (fo) of the 

second frequency (f|p) due to Doppler-effects, 

defining a Doppler frequency interval (fiF-min - fiF-max) around 

the second frequency (fir), the Doppler frequency interval (fiPHnin - 

fiF-max) having a lowest frequency limit (fiF-mm) equal to the 
15 difference between the second frequency (fir) and the maximum 

frequency variation (fo). and a highest frequency limit (fiF-max) equal 

to the sum of the second frequency (fip) and the maximum 

frequency variation (fo). 

dividing the Doppler frequency interval (f|F.min - fiF-max) into an 
20 integer number of equidistant (AO frequency steps (fiF-mfn. fiF-min'*'Af . 

— I fiF-max). and 

defining a frequency candidate vector (fiF-c) for each fre- 
quency step (fiF-fliln. fjF.mln"»*Af, ... , fiFwnax)- 

6. A method according to claim 5, characterized by the steps 
25 of: determining an integer number of initial phase positions 

(9o, 97) for the intermediate frequency signal (Sif). and 

defining a carrier frequencyrphase candidate vector 
(Vf9(fiF.ci ^c)) for each combination of carrier frequency 
candidate vector (fiF.c) and initial phase position (9o. .... Vy). 
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7. A method according to claim 6, characterized by the steps 
of: determining the number of elements (Sn) in each carrier 
frequency-phase candidate vector (yf,(fiF.c. <Pc)). and 
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storing the carrier frequency-phase candidate vector 
(Vf,(ftF-c. ^c)) according to a data format being adapted to per- 
forming correlation operations between the at least one vector 
(Sip.i(k); SiF.Q(k)) derived from the data word (d(k)) and a 
5 segment of a carrier frequency-phase candidate vector (Vf^CfiF-c. 
<Pc)). 

8. A method according to claim 7, characterized by the 
adapting of the data format involving adding at least one ele- 
ment to each segment of the carrier frequency-phase candidate 

10 vector (Vf<p(fiF.c. ^c)) such that the segment attains a number of 
elements which is equal to the number of elements in each the 
at least one vector (SiF.)(k): SiF^(k)) derived from the data word 
(d(k)). thus enabling a segment of the carrier frequency-phase 
candidate vector (Vf<p(fiF.ct ^c)) and one of the at least one 

15 vector (SiF.i(k); StF.Q(k}) to be processed jointly by at least one 
of a SiiyAO-operatioh and an XOR-operation. 

9. A method according to any one of the preceding claims, 
characterized by the steps of: 

determining a maximum variation a code rate due to 
20 Doppler-effects, 

defining a Doppler rate Interval. (CRc-min - CRc.max) around 
a center code rate (CRc). the Doppler frequency interval having 
. a lowest code rate limit (CRcnnin) equal to the difference 
between the center code rate (CRc) and the maximum code rate 
25 . variation (CRq). and a highest frequency limit (CRc-max) equal to 
the sum of the center code rate (CRc) and the maximum code 
rate variation (CRq). 

dividing the Doppler ratei interval (CRcnnin - CRc.max) Into 
an integer number of equidistant (ACR) code rate steps (CRc-min. 

30 CRc-mln+ACR CRc^max)and 

defining a code rate candidate (CRc-c) for each code rate 

step (CRc-min. CRc-min+ACR. .... CRc-max)- 
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10. A method according to claim 9, characterized by deter- 
mining an integer number of possible initial code phase positions 
(0.0 0.9) for each code rate candidate (CRc-c)- 

11. A method according to claim 10, characterized by 
defining, for each signal source specific code sequence (CS(i)). 
a set of combinations between code rate candidate (CRc-c) and 
code phase position (0.0, .... 0.9), each combination repre- 
senting a code rate-phase candidate vector. 

12. A method according to claim 11, characterized by 
generating, for each signal source specific code sequence 
(CS(l)), a set of code vectors (CV) by: 

sampling each code rate-phase candidate vector with the 
basic sampling rate (rs) whereby a corresponding code vector 
(CV) is produced. 

13. A method according to any one of the preceding claims, 
characterized by generating a modified code . vector (CV^) on 
basis of each code vector (CV) by: 

copying a particular number of elements (Ee) from the end 
of an original code vector (CV) to the beginning of the modified 
code vector (CVm), and 

copying the particular number of elements (Eb) from the 
beginning of the original code vector (CV) to the end of the code 
vector (CVm). 

14. A method according to claim 13, characterized by storing, 
for each signal source specific code sequence (CS(i)), a set of 
modified code vectors ({CVm(CRc^, Cph)}), where 

each modified code vector (CV^) contains a number of 
elements (si, .... Sm) representing a sampled version of at least 
one full code sequence (CS). 

a particular modified code vector (CV„,) is defined for each 
combination of code rate candidate (CRc-c) and code phase 




position (Cph). 

15. A method according to ciaim 14, characterized by 
adapting the data format of the modified code vectors (CVm) with 
respect to the data format of the at least one vector (Sif.i(I(); 

5 Sif.q(I()) derived from the data word (d(l<». such that a modified 
code vector (CVm) and one of the at least one vector (SiF.i(ic); 
SiF.Q(i()) may be processed jointly by at least one of a SiMD- 
operation and an XOR-operation. 

16. A method according to any one of the claims 14 or 15, 
10 characterized by involving an initial acquisition phase and a 

subsequent tracking phase, wherein during the acquisition phase 
a set of preliminary parameters are established whi(^ are re- 
quired for initiating a decoding of signals received during the 
tracking phase, a successful acquisition phase resulting in at 
15 least one signal source specific code sequence (CS) being identi- 
fied and a transmitted signal being rendered possible to track, 
each of the at least one signal source specific code sequence 
. (CS) being associated with tracking characteristics in the form of: 
a modified code vector (CVn,), 
20 a carrier frequency candidate vector (fiF-c)> 

an initial phase position (9c), 
a code phase position (Cph), and 

a code index (CI) denoting a starting sample value for the 
modified code vector (CVm). 

25 17. A method according to claim 16, characterized by the 
tracking involving: 

calculating, based the tracking chiaracteristics. a prompt 
pointer (Pp) for each modified code vector (CVm), the prompt 
pointer (Pp) indicating a code sequence start position, and an 
30 initial prompt pointer (Pp) being equal to the code index (CI), and 
assigning, around each prompt pointer (Pp), at least one 
pair of early- and late pointers (Pe, Pl, Pei, Pli Pe2. Plz). where 
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the early pointer (Pe) specifies a sample value being positioned 
at least one element before the prompt pointer's (Pp) position, 
and the late pointer (PO specifies a sample value being 
positioned at least one element after the prompt pointer s (Pp) 
5 position. 

18. A method according to claim 17, characterized by the 

tracking Involving: 

receiving a sequence of Incoming level-discrete sample 

values (1210), . , 

10 forming data words (d(1) d(N)) of the sample values 

(1210) such that each data word (d(k)) contains a number of 
elements which is equal to the number of elements (s„) in each 
carrier frequency-phase candidate vector (Vf,(fiF-ci *Pc)). 

calculating a relevant set of carrier frequency-phase candl- 
15 date vectors (V„(f,F.c. 1»c)) for the data word (d(k)), and 

acquiring, for each carrier frequency-phase candidate 
vector (Vf,(f,F-c, ^c)) in the relevant set. a pre-generated in- 
phase representation (fiFi) and a quadrature-phase represen- 
tation (f,FQ) of the vector (Vf,(f,F-c. 'Pc)) respectively. 

20 19. A method according to claim 18. characterized by the 

tracking Involving: 

multiplying each data word (d(k)) with the in-phase repre- 
sentation (f,Fi) of the carrier frequency-phase candidate vector 
(Vf,(fiF-c, <Pc)) in the relevant set to produce a first intermediate- 
25 frequency-reduced information word (SiF-i(k)). 

multiplying each data word (d(k)) with the quadrature- 
phase representation (f.Fa) of the carrier frequency-phase candi- 
date vector (Vf,(f,F-c, ^c)) in the relevant set to produce a second 
intermediate-frequency-reduced Information word (SiF-aCt)). 

30 20. A method according to claim 19, characterized by the 

tracking involving: 

multiplying the first Intermediate-frequency-reduced infor- 
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matlon word (S„.,(k)) with a modified code f 

starting at a position indicated by the prompt pointer (Pp) to 

produce a first prompt-despread symbol string (A,p(k)). 

multiplying the first intermediate-frequency-reduced m or- 

5 matlon word (S.M(k)) with a modified <^<>^Vl^^^nter^^^^^^ 
starting at a position indicated by an early pomter (Pe) to 
produce a first early-despread symbol string (A,E(k)). 

multiplying the first intermediate-frequency-reduced mfor- 

mation word (S„.,(k)) with a modified <^<>d ^T^l orod'S^e 
10 starting at a position indicated by a late pomter (PO to produce 
a first late-despread symbol string (AiUk)), 

multiplying the second intermedlate-frequency-reduced m- 
formatlon word (S,p.Q(k)) with a modified code vector {CV^.p(k)) 
starting at a position Indicated by the prompt pomter (Pp) to 
15 produce a second prdmpt-despread symbol string (AqpCIc)). 

multiplying the second intermedlate-frequency-reduced In- 
formation word (S,p.Q(k)) with a modified code vector (Cy..E(k)) 
starting at a position indicated by the early pointer (Pe) to 
produce a second early-despread symbol string (AoeCk)), ano 
20 multiplying the second Intermedlate-frequency-reduced in- 

formation word (S,p.Q(k)) with a modified code vector (Cym-L(k)) 
starting at a position Indicated by the late pointer (Pl) to 
produce a second late-despread symbol string (AQL(k)). 

21 A method according to claim 20. characterised by the 
25 tracking involving deriving, for each despread symbol string, a 
resulting data word (DR.,p(k). DR.,E(k). DR,,L(k). DR.Qp(k), Dr^e^k), 

. DR-QL(k)). 

22. A method according to claim 21. characterized by deriving 
the resulting data words (DR.ip(k). DR.,E(k). DR^Uk). DR.Qp(k). Dr- 
30 QE(k); DR<»L(k)) by looking up a respective pre-generated value in 
a table (1310). 
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23. A method according to any one of the claims 19 > 22, 
characterized by performing the multiplication between the data 
word (d(k)) and the in-phase representation (fin) of the carrier 
frequency-phase candidate vector (Vf,(f|i:.c. *Pc)) respective 

5 between the data word (d(lc)) and the a quadrature-phase repre- 
sentation (fiFo) of the carrier frequency-phase candidate vector 
(Vf<p(fiF.c. ^c)) by means of at least one of a SIMD-operation and 
an XOR-'operation 

24. A method according to aiiy one of the claims 19 > 23, 
10 . characterized by performing the multiplication between the Inter- 
mediate-frequency-reduced Information words <SiM(k); SiF.Q(k)) 
and the modified code vectors (CV„,-p, CV„.e; CVm-t) by means of 
at least one of a SIMD-operation and an XOR-operation. 

25. A method according to any one of the claims 19 - 24, 
15 characterized by propagating, in connection with completing 

the processing of a current data word (d(ic}) and initiating the 

processing of a subsequent data word (d(k-'-l)): 

a pointer (Pd) indicating a first, sample value of the 

subsequent data word (dCk-t-l)), 
20 a group of parameters describing the relevant set of carrier 

frequency-phase candidate vectors (Vfq,(f|p.c. ^c)}. 
the relevant set of code vectors (CVm), and 
prompt-, early-, and late pointers (Pp, Pe. Pl}- 

26. A computer program directly loadable into the internal 
25 memory of a computer, comprising software for controlling the 

steps of any of the claims 1-25 when said program is run on 
the computer. 

27. A computer readable medium, having a program recorded 
thereon, where the program is to make a computer control the 

30 steps of any of the claims 1 - 25. 
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28. A signal receiver (1500) for receiving navigation data 
signals transmitted in a navigation satellite system comprising: 

a radio front end unit (1510) adapted to receive a conti- 
nuous radio signal (S„f) and in response thereto produce a 
corresponding electrical signal (S.f) comparatively high 

frequency^^^^^^^ unit (1520) adapted to receive the electrical 
signal (S,f) and in response thereto produce a sequence of 
sample values being divided Into data words (d(k)). and 

a digital processor unit (1530) adapted to receive the data 
words (d(k)) and in response thereto demodulate a data signal, 
and including a memory means (1535). Qharacterlzed in that 
a computer program according to claim 26 is loaded into the 
memory means (1535). 



41 



Abstract 

The present Invention relates to processing of spread spectrum 
signals, where a continuous signal of a comparatively high 
frequency is received (1610). This signal Is sampled at a basic 

5 sampling rate whereby a resulting sequence of time discrete 
signal samples is produced, which are in turn quantised into a 
corresponding level-discrete sample value (1620). A plurality of 
data words are formed, which each includes one or more 
consecutive sample values (1630). Information obtained from 

10 these data words is correlated (1640) with at least one 
representation of a signal source specific code sequence, which 
has been pre-generated in the form of a code vector (1600). The 
correlation step specifically involves correlating at least each 
vector in a sub-group of the code vectors with at least one 

15 vector that has been derived from the data word. Thereby 
resulting data Is produced (1650). 

(Fig. 16) 
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